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Abstract. Turbulence is the most common state of astrophysical flows. In typical as- 
trophysical fluids, turbulence is accompanied by strong magnetic fields, which has a 
large impact on the dynamics of the turbulent cascade. Recently, there has been a signif- 
icant breakthrough on the theory of magnetohydrodynamic (MHD) turbulence. For the 
first time we have a scaling model that is supported by both observations and numerical 
simulations. We review recent progress in studies of both incompressible and compress- 
ible turbulence. We compare froshnikov-Kraichnan and Goldreich-Sridhar models, and 
discuss scalings of Alfven, slow, and fast waves. We also discuss the completely new 
regime of MHD turbulence that happens below the scale at which hydrodynamic tur- 
bulent motions are damped by viscosity. In the case of the partially ionized diffuse 
interstellar gas the viscosity is due to neutrals and truncates the turbulent cascade at 
~parsec scales. We show that below this scale magnetic fluctuations with a shallow 
spectrum persist and discuss the possibility of a resumption of the MHD cascade after 
ions and neutrals decouple. We discuss the implications of this new insight into MHD 
turbulence for cosmic ray transport, grain dynamics, etc., and how to test theoretical 
predictions against observations. 



1 Introduction 

Most astrophysical systems, e.g. accretion disks, stellar winds, the interstellar 
medium (ISM) and intercluster medium are turbulent with an embedded mag- 
netic field that influences almost all of their properties. This turbulence which 
spans from km to many kpc (see discussion in p,15C,36|) holds the key to many 



astrophysical processes (e.g., transport of mass and angular momentum, star 
formation, fragmentation of molecular clouds, heat and cosmic ray transport, 
magnetic reconnection) . Statistics of turbulence is also essential for the cosmic 
microwave background (CMB) radiation foreground studies |8S[ ]. 

All turbulent systems have one thing in common: they have a large "Reynolds 
number" (Re = LV/v\ L= the characteristic scale or driving scale of the sys- 
tem, V=the velocity difference over this scale, and i/=viscosity) , the ratio of the 
viscous drag time on the largest scales (L 2 jv) to the eddy turnover time of a 
parcel of gas (L/V). A similar parameter, the "magnetic Reynolds number" , Rm 
(= LV/r]; ?7=magnetic diffusion), is the ratio of the magnetic field decay time 
(L 2 /rj) to the eddy turnover time (L/V). The properties of the flows on all scales 
depend on Re and Rm. Flows with Re < 100 are laminar; chaotic structures 
develop gradually as Re increases, and those with Re ~ 10 3 are appreciably less 
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chaotic than those with Re ~ 10 7 . Observed features such as star forming clouds 
and accretion disks are very chaotic with Re > 10 8 and Rrn > 10 16 . 

Let us start by considering incompressible hydrodynamic turbulence, which 
can be described by the Kolmogorov theory Suppose that we excite fluid 
motions at a scale L. We call this scale the energy injection scale or the largest 
energy containing eddy scale. For instance, an obstacle in a flow excites motions 
on scales of the order of its size. Then the energy injected at the scale L cascades 
to progressively smaller and smaller scales at the eddy turnover rate, i.e. t» ~ 
vi/l, with negligible energy losses along the cascade []. Ultimately, the energy 
reaches the molecular dissipation scale Id, i.e. the scale where the local Re ~ 1, 
and is dissipated there. The scales between L and Id are called the inertial range 
and it typically covers many decades. The motions over the inertial range are 
self-similar and this provides tremendous advantages for theoretical description. 

The beauty of the Kolmogorov theory is that it does provide a simple scaling 
for hydrodynamic motions. If the velocity at a scale I from the inertial range 
is vi, the Kolmogorov theory states that the kinetic energy (pvf ~ vf as the 
density is constant) is transferred to next scale within one eddy turnover time 
(l/vf). Thus within the Kolmogorov theory the energy transfer rate (vf /(l/vi)) 
is scale-independent, 



constant, (1) 



teas (l/vi) 

and we get the famous Kolmogorov scaling 

vi oc l 1 ^. (2) 

The one-dimensionalQ energy spectrum E(k) is the amount of energy between 
the wavenumber k and k+dk divided by dk. When E(k) is a power law, kE(k) is 
the energy near the wavenumber k ~ l/l. Since vf ~ kE(k), Kolmogorov scaling 
implies 

E(k) oc k~ 5/3 . (3) 

Kolmogorov scalings were the first major advance in the theory of incompress- 
ible turbulence. They have led to numerous applications in different branches of 
science (see | {L12| ). However, astrophysical fluids are magnetized and the a dy- 
namically important magnetic field should interfere with eddy motions. 

Paradoxically, astrophysical measurements are consistent with Kolmogorov 
spectra (see LPE02 p6[]). For instance, interstellar scintillation observations in- 
dicate an electron density spectrum is very close to —5/3 for 10 8 cm - 10 15 cra 
(see Q). At larger scales LPE02 summarizes the evidence of —5/3 velocity power 
spectrum over pc-scales in HI. Solar-wind observations provide in-situ measure- 
ments of the power spectrum of magnetic fluctuations and Leamon et al. [^3| also 
obtained a slope of « —5/3. Is this a coincidence? What properties is the mag- 
netized compressible ISM expected to have? We will deal with these questions, 
and some related issues, below. 

1 This is easy to see as the motions at the scales of large eddies have Re S> 1. 

2 Dealing with observational data, e.g. in LPE02 Js6| , we deal with three dimensional 
energy spectrum P(k), which, for isotropic turbulence, is given by E(k) = ink 2 P(k). 
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Our approach here is complementary to that in Vazquez-Semadeni (this vol- 
ume) and Mac Low (this volume). These reviews discuss attempts to simulate 
the turbulent ISM in all its complexity by including many physical processes 
(e.g. heating, cooling, self-gravity) simultaneously. This provides a possibility of 
comparing observations and simulations (see review by Ostriker, this volume). 
The disadvantage is that such simulations cannot distinguish between the con- 
sequences of different processes. Note, that in studies of turbulence the adaptive 
mesh does not help as the fine structures emerge through the entire computa- 
tional volume. 

Here we discuss a focused approach which aims at obtaining a clear under- 
standing on the fundamental level, and considering physically relevant compli- 
cations later. The creative synthesis of both approaches is the way, we think, 
that studies of astrophysical turbulence should proceed^. Certainly an under- 
standing of MHD turbulence in the most ideal terms is a necessary precursor to 
understanding the complications posed by more realistic physics and numerical 
effects. For review of general properties of MHD, see a recent book by Biskamp 

i- 

In what follows, we first consider observational data that motivate our study 
(§2), then discuss theoretical approaches to incompressible MHD turbulence (§3). 
In §4 we discuss testing and extending of the Goldreich-Sridhar theory of turbu- 
lence, then in §5 we deal with viscous damping of incompressible turbulence and 
describe a new regime of MHD turbulence that is present below the viscous cut- 
off scale. We move to the effects of compressibility in §6 and discuss implications 
of our new understanding of MHD turbulence for the problems of dust motion, 
cosmic ray dynamics, support of molecular clouds, heating of ISM etc in §7. We 
propose observational testing of our results in §8 and present the summary in 
§9. 



2 Observational Data 

Kolmogorov turbulence is the simplest possible model of turbulence. Since it 
is incompressible and not magnetized, it is completely specified by its velocity 
spectrum. If a passive scalar field, like "dye particles" or temperature inhomo- 
geneities, is subjected to Kolmogorov turbulence, the resulting spectrum of the 



passive scalar density is also Kolmogorov (see |95jJl7IJ). In compressible and 
magnetized turbulence this is no longer true, and a complete characterization of 
the turbulence requires not only a study of the velocity statistics but also the 
statistics of density and magnetic fluctuations. 

3 Potentially our approach leads to an understanding of the relationship between mo- 
tions at a given time at small scales (subgrid scales) and the state of the flow at a 
previous time at some larger, resolved, scale. This could lead to a parametrization 
of the subgrid scales and to large eddy simulations of MHD. 
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Direct studies of turbulence^] have been done mostly for interstellar medium 
and for the Solar wind. While for the Solar wind in-situ measurements are possi- 
ble, studies of interstellar turbulence require inverse techniques to interpret the 
observational data. 

Attempts to study interstellar turbulence with statistical tools date as far 
back as the 



1950s [^9 jG4 ; 1 1 7 , 172 1 and various directions of research achieved 
various degree of success (see reviews by J^ , pl] p| j7^ , [79| , p6[ ) . 



2.1 Solar wind 

Solar wind (see review Q) studies allow pointwise statistics to be measured 
directly using spacecrafts. These studies are the closest counterpart of laboratory 
measurements . 

The solar wind flows nearly radially away from the Sun, at up to about 700 
km/s. This is much faster than both spacecraft motions and the Alfven speed. 
Therefore, the turbulence is "frozen" and the fluctuations at frequency / are 
directly related to fluctuations at the scale k in the direction of the wind, as 
k = 2irf/v, where v is the solar wind velocity |57| ]. 

Usually two types of solar wind are distinguished, one being the fast wind 
which originates in coronal holes, and the slower bursty wind. Both of them show, 
however, /~ 5 ^ 3 scaling on small scales. The turbulence is strongly anisotropic 
(see ^6|) with the ratio of power in motions perpendicular to the magnetic 
field to those parallel to the magnetic field being around 30. The intermittency 
of the solar wind turbulence is very similar to the intermittency observed in 
hydrodynamic flows [Bsl. 



2.2 Electron density statistics 

Studies of turbulence statistics of ionized media (see |157| ) have provided infor- 
mation on the statistics of plasma density at scales 10 B -10 15 cm. This was based 
on a clear unde rstan ding of processes of scintillations and scattering achieved by 
theorists^ (see [ 121 ,49|). A peculiar feature of the measured spectrum (see 



is the absence of the slope change at the scale at which the viscosity by neutrals 
becomes important. 

Scintillation measurements are the most reliable data in the "big power law" 
plot in Armstrong et al. (|] . However there are intrinsic limitations to the scintil- 
lations technique due to the limited number of sampling directions, its relevance 
only to ionized gas at extremely small scales, and the impossibility of getting ve- 
locity (the most important!) statistics directly. Therefore with the data one faces 
the problem of distinguishing actual turbulence from static density structures. 

4 Indirect studies include the line-velocity relationships jrHJ where the integrated ve- 
locity profiles are interpreted as the consequence of turbulence. Such studies do not 
provide the statistics of turbulence and their interpretation is very model dependent. 

5 In fact, the theory of scintillations was developed first for the atmospheric applica- 
tions. 
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Moreover, the scintillation data does not provide the index of turbulence di- 
rectly, but only shows that the data are consistent with Kolmogorov turbulence. 
Whether the (3D) index can be -4 instead of -11/3 is still a subject of intense 
debate ]56 , |121 |. In physical terms the former corresponds to the superposition 
of random shocks rather than eddies. 

Additional information on the electron density is contained in the Faraday 
rotation measures of extragalactic radio sources (see [154 155 1 ) . However, there 
is so far no reliable way to disentangle contributions of the magnetic field and the 
density to the signal. We feel that those measurements may give us the magnetic 
field statistics when we know the statistics of electron density better. 



2.3 Velocity and density statistics from spectral lines 

Spectral line data cubes are unique sources of information on interstellar tur- 
bulence. Doppler shifts due to supersonic motions contain information on the 
turbulent velocity field which is otherwise difficult to obtain. Moreover, the sta- 
tistical samples are extremely rich and not limited to discrete directions. In 
addition, line emission allows us to study turbulence at large scales, comparable 
to the scales of star formation and energy injection. 

However, the problem of separating velocity and density fluctuations within 
HI data cubes is far from trivial j7^,(79]]84|,|86| . The analytical description of the 
emissivity statistics of channel maps (velocity slices) in Lazarian & Pogosyan 
[ |S4"| (see also (7^^] for reviews) shows that the relative contribution of the 
density and velocity fluctuations depends on the thickness of the velocity slice. 
In particular, the power-law asymptote of the emissivity fluctuations changes 
when the dispersion of the velocity at the scale under study becomes of the order 
of the velocity slice thickness (the integrated width of the channel map). These 
results are the foundation of the Velocity-Channel Analysis (VCA) technique 
which provides velocity and density statistics using spectral line data cubes. The 
VCA has been successfully tested using data cubes obtained via compressible 
magnetohydrodynamic simulations and has been applied to Galactic and Small 
Magellanic Cloud atomic hydrogen (HI) data p7p4 ,15S,28[. Furthermore, the 



inclusion of absorption effects |85| has increased the power of this technique. 
Finally, the VCA can be applied to different species (CO, H a etc.) which should 
further increase its utility in the future. 

Within the present discussion a number of results obtained with the VCA 
are important. First of all, the Small Magellanic Cloud (SMC) HI data exhibit a 
Kolmogorov-type spectrum for velocity and HI density from the smallest resolv- 
able scale of 40 pc to the scale of the SMC itself, i.e. 4 kpc. Similar conclusions 
can be inferred from the Galactic data [jsD for scales of dozens of parsecs, al- 
though the analysis has not been done systematically. Deshpande et al. |^8| 
studied absorption of HI on small scales toward Cas A and Cygnus A. Within 
the VCA their results can be interpreted as implying that on scales less than 
1 pc the HI velocity is suppressed by ambipolar drag and the spectrum of density 
fluctuations is shallow P(k) ~ fc~ 2,8 . Such a spectrum j27j can account for the 
small scale structure of HI observed in absorption. 
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2.4 Magnetic field statistics 

Magnetic field statistics are the most poorly constrained aspect of ISM turbu- 
lence. The polarization of starlight and of the Far-Infrared Radiation (FIR) from 
aligned dust grains is affected by the ambient magnetic fields. Assuming that 
dust grains are always aligned with their longer axes perpendicular to magnetic 
field (see the review [j8(J), one gets the 2D distribution of the magnetic field 
directions in the sky. Note that the alignment is a highly non-linear process 
in terms of the magnetic field and therefore the magnetic field strength is not 
available^. 

The statistics of starlight polarization (see jio]]) is rather rich for the Galactic 
plane and it allows to establish the spectrum^ E{K) ~ K^ 1 - 5 , where if is a two 
dimensional wave vector describing the fluctuations over sky patchj^j 

For uniformly sampled turbulence it follows from Lazarian & Shutenkov |83[ 
that E(K) ~ K a for K < K and R- 1 for K > K , where K^ 1 is the critical 
angular size of fluctuations which is proportional to the ratio of the injection 
energy scale to the size of the turbulent system along the line of sight. For 
Kolmogorov turbulence a = —11/3. 

However, the real observations do not uniformly sample turbulence. Many 
more close stars are present compared to the distant ones. Thus the interme- 
diate slops are expected. Indeed, Cho & Lazarian |22) showed through direct 



simulations that the slope obtained in 1 40 is compatible with the underlying 



Kolmogorov turbulence. At the moment FIR polarimetry does not provide maps 
that are really suitable to study turbulence statistics. This should change soon 
when polarimetry becomes possible using the airborne SOFIA observatory. A 
better understanding of grain alignment (see pQ|) is required to interpret the 
molecular cloud magnetic data where some of the dust is known not to be aligned 
(see |$2j and references therein). 

Another way to get magnetic field statistics is to use synchrotron emission. 
Both polarization and intensity data can be used. The angular correlation of po- 
larization data H shows the power-law spectrum K~ 18 and we believe that the 
interpretation of it is similar to that of starlight polarization. Indeed, Faraday 
depolarization limits the depth of the sampled region. The intensity fluctuations 
were studied in with rather poor initial data and the results were inconclu- 
sive. Cho & Lazarian [p2| interpreted the fluctuations of synchrotron emissivity 
[EMB in terms of turbulence with Kolmogorov spectrum. 



The exception to this may be the alignment of small grains which can be revealed 
by microwave and UV polarimetry |8C(| . 

7 Earlier papers dealt with much poorer samples (see [[i5|) and they did not reveal 
power-law spectra. 

8 This spectrum is obtained by j!o| in terms of the expansion over the spherical 
harmonic basis Y; m . For sufficiently small areas of the sky analyzed the multipole 
analysis results coincide with the Fourier analysis. 
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3 Theoretical Approaches to MHD Turbulence 



Here we consider mainly Kolmogorov-type theories. Other theories not discussed 
in this sectio n inc lude the eddy-damped quasinormal Markovian (EDQNM) ap- 
proximation [139 1, the renormalization group technique |3£, 174 |, and the direct 
interaction approximation f70fl. 



3.1 Iroshnikov-Kraichnan theory 

Attempts to describe magnetic turbulence statistics were made by Iroshnikov 
[ |6l| and Kraichnan JnJ . Their model of turbulence (IK theory) is isotropic in 
spite of the presence of the magnetic field. 

We can understand the IK theory as follows .[] For simplicity, let us suppose 
that a uniform external magnetic field (Bo) is present. In the incompressible 
limit, any magnetic perturbation propagates along the magnetic field line. To 
the first order, the speed of propagation is constant and equal to the Alfven 
speed Va = Bo/^Anp, where p is the density. Since wave packets are moving 
along the magnetic field line, there are two possible directions for propagation. 
If all the wave pac kets are moving in one direction, then they are stable to 



nonlinear order 1 1 32 1. Therefore, in order to initiate turbulence, there must be 
opposite-traveling wave packets and the energy cascade occurs only when they 
collide. The IK theory starts from this observation. 

The IK theory assumes that, when two opposite-traveling wave packets of 
size I collide, they lose the following amount of energy to smaller scales: 

AE ~ {dv 2 /dt)At ~ v, ■ viAt ~ vi{vf/l)At ~ (vf/l)(l/V A ), (4) 

where the IK theory assumes that only collisions between similar size packets are 
important and At ~ I/Va- The latter means that the duration of the collision 
is the size of the wave packet divided by the speed of the wave packets. These 
assumptions look reasonable at first. But, it is important to note that they fail 
when eddies are anisotropic. That is, if eddies are elongated along the magnetic 
field line, then At is not I/Va, but 1\\/Va, where l\\ is the parallel size of the wave 
packet (or 'eddy'). 

Equation (f|) tells us that the energy change per collision is vf (vi/Va), which 
is only a tiny fraction of vf when Va 3> «z- Therefore, in order for the eddy to 
transfer all the energy to small eddies, the eddy must go through many collisions. 
When such collisions are incoherent, we require a total (v 2 /AE) 2 collisions to 
complete the cascade. This means that the energy cascade time t cas is 

\ AE J vi vi 

which means that this new cascade time is (Va/vi) times longer than the eddy 
turnover time (l/vi). As in the Kolmogorov theory, the IK theory assumes the 

9 We follow arguments in ]46[. 
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constancy of energy cascade (eq. : {vf)/ {IVa) — constant, which, in turn, yields 
vi oc Z 1 / 4 , or, 

Iroshnikov-Kraichnan: E{k) oc fc~ 3 / 2 . (6) 

3.2 Anisotropy 

A uniform component to the magnetic field defines a special direction, which will 
be reflected in the dynamics of turbulent fluctuations. One obvious effect is that 
it is easy to mix field lines in directions perpendicular to the local mean magnetic 
field and much more difficult to bend them. The IK theory assumes isotropy of 
the energ y cascade in Fourier space, an assumption which has attracted severe 



criticism [115,152,114,158 104 . Mathematically, anisotropy manifests itself in 



the resonant conditions for 3-wave interactions: 

k 1 +k a =k 3 , (7) 

U\ + ui 2 = ^3, (8) 

where k's are wavevectors and w's are wave frequencies. The first condition is 
a statement of wave momentum conservation and the second is a statement of 
energy conservation. Alfven waves satisfy the dispersion relation: u> = y^fey], 
where fcii is the component of wavevector parallel to the background magnetic 
field. Since only opposite-traveling wave packets interact, ki and k2 must have 
opposite signs. Then from equations (Q) and (ft), either fen x or fen 2 must be 
equal to and fen 3 must be equal to the nonzero initial parallel wavenumber. 



That is, zero frequency modes are essential for energy transfer [153]. Therefore, 
in the wavevector space, 3-wave interactions produce an energy cascade which is 
strictly perpendicular to the mean magnetic field. However, in real turbulence, 
equation (||) does not need to be satisfied exactly, but only to within an an error 
of order 5uj ~ l/t cas pcf . This implies that the energy cascade is not strictly 
perpendicular to Bo, although clearly very anisotropic. 

It is noteworthy that there has been a claim that "by increasing the mag- 
nitude of the mean field in 3-D simulations one finds that the transition from 
th e iso tropic 3-D scaling properties toward those observed in 2-D" . (see also 
[([ 175, (| for recent development in 2-D MHD turbulence). This claim has yet 



to be substantiated, however. We feel that the available numerical simulations 
[ 116 2^,101 1 are reasonably consistent with the Goldreich & Sridhar |^5| model 



that we review in the next section. It is also worth noting that the idea of an 
anisotropic (perpendicular) cascade has been incorporated into the framework 
of the reduced MHD approximation ]i6f|,|f46|,pl3|Jl8^fl . 



3.3 Goldreich-Sridhar theory 

We assume throughout this discussion that the rms turbulent velocity at the 
energy injection scale is comparable to the Alfven speed of the mean field and 
consider only scales below the energy injection scale. Consequently, we are not 
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concerned with the problem of magnetic field generation, or the magnetic dy- 
namo which is considered elsewhere in this volume 

(see also [111 132,72]] for reviews, 107 , 123| , 12, 15 for numerical calculations, 



[ p73[p[l5t for suppression of the a dynamo effect in highly conducting fluids 



and [73,11 176 100 , 185 for recent developments). 

An ingenious model very similar in its beauty and simplicity to the Kol- 
mogorov model has been proposed by Goldreich & Sridhar [[l5| (1995; hereinafter 
GS95) for incompressible MHD turbulence. They pointed out that motions per- 
pendicular to the magnetic field lines mix them on a hydrodynamic time scale, 
i.e. at a rate t~^ s ss k±vi, where k± is the wavevector component perpendicular 
to the local mean magnetic field and I ~ /c _1 (ss fcj 1 )- These mixing motions cou- 
ple to the wave-like motions parallel to magnetic field giving a critical balance 
condition 

k\\V A ~ k±v k , (9) 

where fcii is the component of the wavevector parallel to the local magnetic field. 
When the typical fc|| on a scale k± falls below this limit, the magnetic field tension 
is too weak to affect the dynamics and the turbulence evolves hydrodynamically, 
in the direction of increasing isotropy in phase space. This quickly raises the 
value of fen. In the opposite limit, when fc|| is large, the magnetic field tension 
dominates, the error Sui in the matching conditions is reduced, and the nonlinear 
cascade is largely in the k± direction, which restores the critical balance. 

If conservation of energy in the turbulent cascade applies locally in phase 
space then the energy cascade rate (vf /t cas ) is constant (eq. |l|): (vf)/(l/vi) = 
constant. Combining this with the critical balance condition we obtain an 
anisotropy that increases with decreasing scale 

cx k 2 / 3 , (10) 

and a Kolmogorov-like spectrum for perpendicular motions 

vi cx Z 1/3 , or, E(k) cx fc^ 5/3 , (11) 

which is not surprising since the magnetic field does not influence motions that 
do not bend it. At the same time, the scale-dependent anisotropy reflects the 
fact that it is more difficult for the weaker, smaller eddies to bend the magnetic 
field. 

GS95 shows the duality of motions in MHD turbulence. Those perpendicular 
to the mean magnetic field are essentially eddies, while those parallel to magnetic 
field are waves. The critical balance condition couples these two types of motions. 

3.4 Weak/Intermediate turbulence. 

Let us reconsider the interaction of two wave packets moving oppositely along 
the mean magnetic field line. As in equation (Q), the energy loss per collision is 

AE ~ (v*)At ~ (vfAOGy/VA), (12) 
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where we explicitly distinguish the parallel (7||) and the perpendicular size (lx)- 
The kinetic energy of the eddy is vf. Therefore, the ratio of AE to E, 

r = AE v h _ vik± n ^ 
° " vf ~ V A l x ~ V A k {l ' 

characterizes the strength of the nonlinear interaction [[l5| . In the GS95 theory, 
~ 1. This means that Va ~ V is required at the energy injection scale when 
energy injection is isotropic (fcj_.L ~ fell, l)- When this condition is satisfied, the 
turbulence is called strong turbulence. 



There are some astrophysical situations, e.g. the Jovian magnetosphere [148], 
where the parameter £j, is much smaller than the unity over a broad range of 
length scales. Although as noted above the cascade will evolve in the direction of 
increasing Q for decreasing I, and may reach the strong turbulent regime on very 
small scales. In this regime, the parallel cascade is strongly suppressed so the 
turbulence is qualitatively different from the strong turbulence discussed above. 
This is the weak turbulence regime. We do not discuss this type of turbulence here 
due to its limited astrophysical applicability and restricted inertial range. For 
more information, see |l^ , |l58|Jl22p ^,^| . Note that Galtier et al. Q obtained 
E(k) ~ kj 2 (see also |46| 



4 Testing and Extending Incompressible Theory 

Here we focus on recent direct numerical simulations related to the anisotropic 
structure of MHD turbulence. A discussion of earlier pioneering numerical sim- 



ulations of MHD turbulence can be found in [140 107 1. 



4.1 Scaling laws 

Despite its attractiveness, the Goldreich-Sridhar model is a conjecture that re- 
quires testing. The first such test was done by Cho & Vishniac Q who used an 
incompressible pseudo-spectral MHD code with 256 3 grid points. In their simu- 
lations, they used Va — B / y/Anp ~ V. Their results for eddy shapes are shown 
in Figure Qa, which shows a reasonable agreement with the predicted scale- 
dependent anisotropy of the turbulence (eq. (|Io|)) in the inertial range (i.e. the 
scales between the energy injection scale and the dissipation scale). They also 
obtained E(k) ~ fc~ 5 / 3 in the inertial range. Yet although the velocity does show 
the expected scaling, the magnetic field scaling is a bit more uncertain. 



A subsequent numerical study by Maron & Goldreich [ 1 1 1 performed with 
a different code and in a different physical regime, namely, for Va ^> V, also 
supported the GS95 model and clarified the role of pseudo-Alfven and Alfven 
modes. In particular, they confirmed that the pseudo-Alfven modes are passively 
carried down the cascade through interactions with the Alfven modes. They also 
showed that passive scalars adopt the same power spectrum as the velocity and 
magnetic field fluctuations. In addition, they addressed several issues about the 
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Fig. 1. (a) Eddies and local mean magnetic field. Local mean field is the properly 
averaged magnetic field near eddies under consideration. Eddies at different locations 
(e.g. eddies 1 & 1') can have different local mean fields. Eddies of different sizes (e.g. ed- 
dies 1 & 2) can also have different local mean fields. When we do not take into account 
the local mean field, calculations hardly reveal true eddy structures. From ps] ]. (b) An 
example of calculating the local mean field. The direction of the local mean field B( is 
obtained by the pair-wise average: B; = (B(ri) + B(r2))/2. 



imbalanced cascade. Overall, they obtained results reasonably consistent with 
the GS95 model, but their energy spectra scale slightly differently: E(k) oc fc~ 3 / 2 . 
They attributed this result to intermittency. 

Mullcr & Biskamp [116] studied MHD turbulence numerically in the regime of 
V 3> Va and obtained a Kolmogorov spectrum: E{k) ~ fc~ 5 / 3 . They numerically 
studied the scaling exponents and obtained £2 ~ 0.7 and C3 ~ 1 (see the definition 
of the scaling exponents in the next subsection). The value of £2 is consistent 
with the energy spectrum. 

Other related rece nt numerical simu latio ns include Matthaeus et al. [104 
and Milano et al. [10£|. Matthaeus et al. [104] showed that the anisotropy of low 
frequency MHD turbulence scales linearly with the ratio of perturbed and total 
magnetic field strength b/B (= b/(b 2 + £ 2 ) 1/2 ).0 



In fact we can derive the GS95 scaling using this result. Although their analysis 
was based on comparing the strength of a uniform background field and the mag- 
netic perturbations on all scales, we can reinterpret this result by assuming that 
the strength of random magnetic field at a scale I is bi, and that the background 
field is the sum of all contributions from larger scales. Then Matthaeus et al.'s re- 
sult becomes a prediction that the anisotropy (fc||/fcj_) is proportional to (bi/B). We 
can take the total magnetic field strength B ~ constant as long as the background 
field is stronger than the perturbations on all scales. Since bi ~ (kE(k)) 1 ^ 2 ~ fc^ 1//3 , 
we obtain an anisotropy (k\\/k±) proportional to k^^ 3 , and k\\ oc k 2 / 3 . In this in- 
terpretation, smaller eddies are more elongated because they have a smaller bi/B 
ratio. 




Fig. 2. Cross-sections of the data cube. (Left panels) |b;| in a plane || to Bq. (Right 
panels) |b; | in a plane _L to Bo. In the left panels, Bo is along the horizontal axis. Large 
scale eddies are obtained from the Fourier components with 1 < k < 4. Medium scale 
eddies are obtained from the Fourier components with 4 < k < 16. Small scale eddies 
are obtained from the Fourier components with 16 < k < 64. The small scale eddies 
show a high degree of elongation in the parallel plane. However, they do not show a 
systematic behavior in the perpendicular plane. 



All in all, numerical simulations so far have been largely, but not perfectly, 
consistent with the GS95 theory, e.g. the Kolmogorov-type scaling and the scale- 
dependent anisotropy (fc|| oc k 2 / 3 ), and helped to extend it. An important point, 
which was not included in some of the earlier work, is that the scale dependent 
anisotropy can be measured only in a local coordinate frame which is aligned 
with the locally averaged magnetic field direction ]2Q| . The necessity of using 
a local frame is due to the fact that eddies are aligned along the local mean 
magnetic fields, rather than the global mean field Bq. Since smaller scale eddies 




Fig. 3. (a) Velocity correlation function (VCF) from a simulation. The contours repre- 
sent shape of different size eddies. The smaller contours (or, eddies) are more elongated, 
(b) VCF generated from equation (Q. From Q. 
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Fig. 4. (a) Semi-major axis and semi-minor axis of contours of the VCF shown in 
Figure 0a. The results support kn oc k 2 / 3 . From @. (b) Scaling exponents. From |[23|] . 



are weaker, and more anisotropic, measurements of eddy shape based on a global 
coordinate system are always dominated by the largest eddies in the simulation. 
Figure [l] illustrates the concept of a local frame and one way to identify it. 
Further research in Cho, Lazarian, & Vishniac [p3| (2002a: hereinafter CLV02a) 
showed that in the local system of reference the mixing motions perpendicular 
to the magnetic field have statistics identical to hydrodynamic turbulence (cf. 
Miiller & Biskamp pf|). 

Fig. U shows the shapes of eddies of different sizes. Left panels show an 
increased anisotropy as we move from the top (large eddies) to the bottom (small 
eddies). The horizontal axes of the left panels are parallel to Bo- Structures in 
the perpendicular plane do not show a systematic elongation. 

Fig. ||a and Fig. [i]a quantify some of these results. The contours of the corre- 
lation function obtained in 2f| are shown in Fig. tm and are consistent with the 
predictions of the GS95 model. Fig. ||a shows that the semi-major axis (l/fc||) 
is proportional to the 2/3 power of the semi-minor axis (1/fc^), implying that 

2/3 

fc|| OC kl . While the one dimensional energy spectrum follows Kolmogorov spec- 
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trum, E(k) oc k 5 / 3 , CLV02a showed that the 3D energy spectrum can be fit 

by 

P(k ± ,k {l ) = (B /L^)kl W/3 e X p (-^ 1/3 ^) , (14) 

where Bq is the strength of the mean field and L is the scale of energy injection. 
The velocity correlation from the 3D spectrum provides an excellent fit to the 
numerical data (Fig. |^b). This allows practical applications illustrated in §5. 



4.2 Intermittency 

Intermittency refers to the non-uniform distribution of dissipative structures. 
Intermittency has an important dynamical consequence: it affects the energy 
spectrum. Highly intermittent turbulent structures were invoked by Falgarone et 
al. ]37j and Joulain et al. |6j| as the primary location of endothermic interstellar 
chemical reactions. 



Maron & Goldreich [101 studied the intermittency of dissipation structures 
in MHD turbulence using the fourth order moments of the Elsasser fields and 
the gradients of the fields. Their simulations show strong intermittent structures. 
CLV02a used a different, but complementary, method to study intermittency, 
based on the higher order longitudinal structure functions. They found that 
by this measure the intermittency of velocity field in MHD turbulence across 
local magnetic field lines is as strong as, but not stronger than, hydrodynamic 
turbulence. 

In fully developed hydrodynamic turbulence, the (longitudinal) velocity struc- 
ture functions S p —< ([v(x + r) — v(x)] ■ f) p >=< 5v p L {r) > are expected 
to scale as p". For example, the classical Kolmogorov phenomenology (K41) 
predicts ( p — p/3. The (exact) result for p=3 is the well-known 4/5-relation: 
< 5v^(j) >= — (4/5)er, where e is the energy injection rate (or, energy dis- 
sipation rate) (see e.g. El)). On the other hand, considering intermittency, 



She & Leveque ([152|; hereinafter S-L) proposed a different scaling relation: 
Cp L = P/9 + 2[1 - (2/3) p/3 ]. Note that the She-Leveque model also implies 

C3 = 1. 

So far in MHD turbulence, to the best of our knowledge, there is no rigor- 
ous intermittency theory which takes into account scale-dependent anisotropy. 
Politano & Pouquet |135[ have developed an MHD version of the She-Leveque 
model: 

*P p = Eh _ ~\ M _ (-\ - rr/nw/s 



Cr = ~(1 - *) + C (1 - (1 - x/Cf'o) , (15) 

where C is the co-dimension of the dissipative structure, g is related to the 
scaling vi ~ l 1 ^ 9 , and x can be interpreted as the exponent of the cascade time 
teas l x • (In fact, g is related to the scaling of Elsasser variable z=v ± b: 
zi ~ I 1 / 9 .) In the framework of the IK theory, where g — 4, x — 1/2, and 
C = 1 when the dissipation structures are sheet-like, their model of intermittency 



becomes C^ K = p/8 + 1 - (l/2) p / 4 . On the other hand, Miiller & Biskamp |TT6| 
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performed numerical simulations on decaying isotropic MHD turbulence and 
obtained Kolmogorov-like scaling (E(k) ~ fc~ 5 / 3 and t ~ I 2 / 3 ) and sheet-like 
dissipation structures, which implies g = 3, C = 1, and x = 2/3. From equation 
(p^|), they proposed that 

C p MS = P /9 + l-(l/3f/ 3 . (16) 

The intermittency results from [^3| are shown in Fig. [|b. The filled circles 
represent the scaling exponents of longitudinal velocity structure functions in 
directions perpendicular to the local mean magnetic field. It is surprising that 
the scaling exponents are so close the original (i.e. hydrodynamic) S-L model. 
This raises an interesting question. In the simulations of CLV02a, t cas oc I 2 ! 3 
and E(k) oc k~ 5 ^ 3 scaling is observed. It is observed that MHD turbulence has 



sheet-like dissipation structures |137. Therefore, the parameters for CLV02a 
simulations should be the same as those of Muller & Biskamp's (i.e. g — 3, C = 1, 
and x — 2/3) rather than suggesting C = 2. We believe this difference stems from 
the different simulation settings (Muller & Biskamp's turbulence is isotropic and 
CLV02a's is anisotropic) and the way the turbulence is analyzed (global versus 
local frame). In fact, we expect the the small scale behavior of MHD turbulence 
should not depend on whether the largest scale fields are uniform or have the 
same scale of organization as the largest turbulent eddies. Nevertheless, given the 
limited dynamical range available in these simulations, it would not be surprising 
if the scale of the magnetic field has an impact on the intermittency statistics. It 
is not clear how scale-dependent anisotropy changes the intermittency model in 
equation ( |l5| ) and we will not discuss this issue further. Instead, we simply stress 
that a striking similarity exists between ordinary hydrodynamic turbulence and 
MHD turbulence in perpendicular directions, which further supports the picture 
of the GS95 turbulence.^ 

In Figure ||a, we also plot the scaling exponents (represented by filled squares) 
of longitudinal velocity structure functions along directions of the local mean 
magnetic field. Although we show only the exponents of longitudinal struc- 
ture functions, those of transverse structure functions follow a similar scaling 
law. Evidently intermittency along the local mean magnetic field directions is 
completely different from the scaling predicted by previous (isotropic) models. 
Roughly speaking, the scaling exponents along the directions of local magnetic 
field arc 1.5 times larger than those of perpendicular directions. This has an 
obvious similarity to the scaling of eddy shapes. 

The second order exponent £2 is related to the the 1-D energy spectra: 
E(k ± ) cx fc^ (1+C2) . Previous 2-D driven MHD calculations for B = by Poli- 



tano, Pouquet, & Carbone |136| also found (2 ~ 0.7. However, Biskamp, & 
Schwarz Q| obtained (2 ~ 0.5 from decaying 2-D MHD calculations with Bq = 0. 
The result of CLV02a suggests that (2 is closer to 2/3, rather than to 1/2. (It 
is not clear whether or not the scaling exponents follow the original S-L model 
exactly. At the same time, our calculation shows that the original S-L model can 

11 MG01 attributed the deviation of their spectrum from the Kolmogorov-type to the 
intermittency present in their simulation. 
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be a good approximation for our scaling exponents. The S-L model predicts that 
(2 ~ 0.696.) This is equivalent to our earlier claim that our result supports the 
scaling law E(k±) cx fc ± 5 ^ 3 at least for velocity. For the parallel directions, the 
results support E(k^) oc k^ 2 although the uncertainty is large. 

5 Damping of Turbulence 
5.1 Imbalanced cascade 

Turbulence plays a critical role in molecular cloud support and star formation 
and the issue of the time scale of turbulent decay is vital for understanding these 
processes. If MHD turbulence decays quickly, then serious problems face the re- 
searchers attempting to explain important observational facts, e.g. turbulent mo- 
tions seen within molecular clouds without star formation [111 91 and rates of star 



formation [ 106 1 . Earlier studies attributed the rapid decay of turbulence to com- 
pressibility effects |)^]. GS95 predicts and numerical simulations, e.g. CLV02a, 
confirm that MHD turbulence decays rapidly even in the incompressible limit. 
This can be understood if mixing motions perpendicular to magnetic field lines 
are considered. As we discussed earlier, such eddies, as in hydrodynamic turbu- 
lence, decay in one eddy turnover time. 



Below we consider the effect of imbalance |103,165, pO[|6C ,101| on the turbu- 
lence decay time scale. Duality of the MHD turbulence means that the turbulence 
can be described by colliding wave packets. 'Imbalance' means that the flux of 
wave packets traveling in one direction is significantly larger than those travel- 
ing in the other direction. In the ISM, many energy sources are localized both 
in space and time. For example, in terms of energy injection, stellar outflows 
are essentially point energy sources. With these localized energy sources, it is 
natural that interstellar turbulence be typically imbalanced. 
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Here we show results of the CLV02a study that demonstrate that imbalance 
does extend the lifetime of MHD turbulence (Fig. ||a). We used a run on a grid 
of 144 3 to investigate the decay time scale. For initial conditions we took a data 
cube from a driven turbulence simulation. The initial data cube contains both 
upward (denoted as +) and downward moving waves (denoted as — ). To adjust 
the degree of initial imbalance, we either increased or decreased the energy of the 
upward moving components components and, by turning off the forcing terms, 
let the turbulence decay. Note that the initial energy is normalized to 1. The 
y-axis is the normalized total (=up + down) energy. 

The dependence of the turbulence decay time on the degree of imbalance 
is an important finding. To what degree the results persist in the presence of 
compressibility is the subject of our current research. It is obvious that results of 
CLV 02a are applicable to incompressible, namely, Alfven motions]^] We show in 



i.3 that the Alfven motions arc essentially decoupled from compressible modes. 



As the result we expect that the turbulence decay time may be substantially 
longer than one eddy turnover time provided that the turbulence is imbalanced. 



5.2 Ion-neutral damping: a new regime of turbulence 

In hydrodynamic turbulence viscosity sets a minimal scale for motion, with an 
exponential suppression of motion on smaller scales. Below the viscous cutoff 
the kinetic energy contained in a wavenumber band is dissipated at that scale, 
instead of being transferred to smaller scales. This means the end of the hydrody- 
namic cascade, but in MHD turbulence this is not the end of magnetic structure 
evolution. For viscosity much larger than resistivity, v ^> 77, there will be a broad 
range of scales where viscosity is important but resistivity is not. On these scales 
magnetic field structures will be created by the shear from non-damped turbu- 
lent motions, which amounts essentially to the shear from the smallest undamped 
scales. The created magnetic structures would evolve through generating small 
scale motions. As a result, we expect a power-law tail in the energy distribution, 
rather than an exponential cutoff. To our best knowledge, this is a completely 
new regime for MHD turbulence. 

In partially ionized gas neutrals produce viscous damping of turbulent mo- 
tions. In the Cold Neutral Medium (see Draine & Lazarian JmJ for a list of the 
idealized phases) this produces damping on the scale of a fraction of a parsec. 
The magnetic diffusion in those circumstances is still negligible and exerts an 
influence only at the much smaller scales, ~ 100km. Therefore, there is a large 
range of scales where the physics of the turbulent cascade is very different from 
the GS95 picture. 

12 In the long run, the imbalance will be defeated by the parametric instability, which 
develops thr ough f ormation of density inhomogeneities within the beam of Alfven 
waves | p^ , ^7|j62[p^ ] . However, this instability takes many wave periods to be estab- 
lished. A similar argument can be applied when we consider completely imbalanced 
cascade. That is, even in the completely imbalanced cascade decay of energy can 
occur due to non-linear steepening if waves. But, this will be very slow for Alfven 
waves. 
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Cho, Lazarian, & Vishniac pij have explored this regime numerically. Here 
we used a grid of 384 3 and a physical viscosity for velocity damping. The kinetic 
Reynolds number is around 100. With this Reynolds number, viscous damping 
occurs around k ~ 7, or about ~ 1/7 of the width of the computational box. 
We minimized magnetic diffusion through the use of a hyper-diffusion term of 
order 3. To test for possible "bottle neck" effects we also did simulations with 
normal magnetic diffusion and reproduced similar results but only over a re- 
duced dynamical range. The bottleneck effect is a common feature in numerical 
hydrodynamic simulations with hyperviscosity. See [^) for MHD simulations. 

In Fig. ||b, we plot energy spectra. The spectra consist of several parts. First, 
the peak of the spectra corresponds to the energy injection scale. Second, for 
2 < k < 7, kinetic and magnetic spectra follow a similar slope. This part is 
more or less a severely truncated inertial range for undamped MHD turbulence. 
Third, the magnetic and kinetic spectra begin to decouple at k ~ 7. Fourth, after 
k ~ 20, a new damped- scale inertial range emerges. In the new inertial range, 
magnetic energy spectrum follows 

New Regime: Eb(k) tx k , (17) 

implying considerable magnetic structures below the viscous damping scale. The 
velocity power spectrum steepens in this regime, but does not fall exponentially. 

Figure ||a shows that the small magnetic structures are highly intermittent 
in the viscous-damped regime. Here we obtained the small scale magnetic field 
by eliminating Fourier modes with k < 20. We can see that the typical radius of 
curvature of field lines in the plane is much larger than the typical perpendicular 
scale for field reversal. The typical radius of curvature of field lines corresponds to 
the viscous damping scale, indicating that stretched structures are results of the 
shearing motions at the viscous damping scale. There is no preferred direction for 
these elongated structures. A similar plot for ordinary (not viscously damped) 
MHD turbulence (Fig. ^b) shows much less intermittent structures. We remove 
Fourier modes with k < 20 also in Fig. ^|b. 

A theoretical model for this new regime and its consequences for stochastic 
reconnection p9| will be found in an upcoming paper (Lazarian, Vishniac, & Cho 
2002). Here we summarize the main points of the model. We begin by noting 
that the strong intermittency seen in this regime suggests a new parameter, 
which is the volume filling fraction of structures with scales comparable to I. This 
in turn implies that we need to distinguish between volume averaged means and 
typical values of velocity and magnetic field perturbations in fraction of space 
where they are concentrated. We denote the latter with a so that vf — fivf, 
and tif = fibf. This model does not include information about the range of field 
strengths or velocities with structure on a scale I in the bulk of the volume, aside 
from assuming that they are sufficiently weak that they do not contribute to any 
global averages. 

The second new fundamental parameter in this model is the eddy turn over 
rate at the damping scale, kd, i.e. 

Tg 1 ~ k d v d ~ k 2 d v n , (18) 
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Fig. 6. {Left) Viscous-damped turbulence. Strength of magnetic field in a plane per- 
pendicular to Bo. Arrows are magnetic fields in the plane. Only a part of the plane is 
shown. Note highly intermittent structures. From 1241] . (Right) Same as Left, but for 
ordinary (not viscous-damped) MHD turbulence. Structures are less intermittent. 



where Vd is the velocity at Id ~ 1/fcd and v n is the viscosity of the plasma due 
to neutral particles, which differs from the viscosity of neutral fluid by the ratio 
of atomic to total densities. 

Since motions on smaller scales are strongly damped, the cascade of magnetic 
energy to smaller scales is due to motions on the damping scale. This implies 
that 

bf/r s ~ constant, (19) 

so that bi is maintained approximately at the level of magnetic field at the 
damping scale, bd- 

This folding and refolding, perpendicular to the mean field direction, has a 
weak effect on the field line curvature. (This result is also seen in the simulations.) 
The decrease in the structure length I is due to an increase in the magnetic 
field gradient perpendicular to the mean field direction. The resulting magnetic 
pressure gradients are balanced by plasma pressure gradients. Thus on the scales 
below the viscous cutoff, the tension forces are balanced by viscous drag, i.e. 

~ xaaxlkjn^Bofa. (20) 

Finally, this dynamic equilibrium can be maintained only if the small scale 
motions are strong enough to counteract the shear, t~ . In other words, vi/l ~ 

Combining these results we see that /; ~ k^l, bi ~ bd(kdl) and vi oc 
/ 3 / 2 . These scaling laws are at least qualitatively consistent with the simulation 
results, although the velocity power spectrum may be slightly steeper than the 
model prediction. 
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Unfortunately, a realistic treatment of the ISM requires an explicit recogni- 
tion of the two fluid nature of the partially ionized plasma, rather than simply 
representing neutral drag with an effective viscosity. Here we are beyond the 
reach of available simulations, and need to rely on an extension of the scaling 
arguments given above. 

First, at sufficiently small values of I, the ambipolar diffusion time (~ 
k 2 U n (b 2 I ptot){Pn/ Pi), where U n = ion-neutral collision time) will become less 
than t Si and the magnetic pressure gradients will be supported entirely by the 
ionized particle pressure. This means that only a fraction, ~ pi/ptot, of the 
energy will continue the cascade to smaller scales. 

Next, for some ion-neutral collision time, ti n , there will be some decoupling 
scale, l c , where 

pVn/ll ~ Pi/tin- (21) 

At smaller scales the ions will be dragged through a more or less uniform neutral 
background. The argument above needs to be modified by replacing the viscous 
drag coefficient in equation ( |20| ) with the right hand side of equation (|2l]). This 
gives 

fi oc r 1 , (22) 

k cx l 1/2 , and vi cx l 1 ' 2 . (23) 

Finally, at some sufficiently small scale the filling factor will rise to unity, 
and the gradients in the magnetic field will become strong enough that neutral 
drag can be ignored. These conditions are satisfied simultaneously when 

k d l~t ln /T s . (24) 

Below this scale we expect to see a resumption of the turbulent cascade, now 
involving only the ionized component of the plasma, down to scales where plasma 
resistivity and viscosity finally dissipate the remaining energy. Since the longest 
cascade time for this regime is £;„ -C t s , we expect this small scale turbulence 
to be intermittent, with a duty cycle ~ U n /T s . 

All the consequences of the new regime of the MHD turbulence have not yet 
been appreciated, but we expect that it will have a substantial impact on our 
understanding of the interstellar physics. Moreover, the treatment given above 
actually applies only when pij p is not very small. Otherwise the decoupling scale 
can be larger than the viscous damping scale. 

6 Compressible Turbulence 

For the rest of the review, we consider MHD turbulence of a single conducting 
fluid. While the GS95 model describes incompressible MHD turbulence well, 
no universally accepted theory exists for compressible MHD turbulence despite 
various attempts (e.g., pq|). Earlier numerical simulations of compressible MHD 
turbulence covered a broad range of astrophysical problems, such as the decay 
of turbulence (e.g. ]9^ , |l6(| ) or turbulent modeling of the ISM (see recent review 
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Fig. 7. (a) Directions of fast and slow basis vectors. £/ and £ 3 represent the directions 
of displacement of fast and slow modes, respectively. In the fast basis (£/) is always 
between k and k±. In the slow basis (£ s ) lies between # and Bo - Here, 9 is perpendicular 
to k and parallel to the wave front. All vectors lie in the same plane formed by Bo 
and k. On the other hand, the displacement vector for Alfven waves (not shown) is 
perpendicular to the plane, (b) Directions of basis vectors for a very small drawn 
in the same plane as in (a). The fast bases are almost parallel to kx. (c) Directions 
of basis vectors for a very high /3. The fast basis vectors are almost parallel to k. The 
slow waves become pseudo- Alfven waves. 



pO) ; see also |3|,|6§(l3|,|l6|l for earlier 2D simulations and |2|,|26|JT3T[^,| 

for recent 3D simulations). In what follows, we concentrate on the fundamental 
properties of compressible MHD. 



6.1 Alfven, slow, and fast modes 

Let us start by reviewing different MHD waves. In particular, we describe the 
Fourier space representation of these waves. The real space representation can 
be found in papers on modern shock-capturing MHD codes (e.g. |l5| , |147| ]). For 
the sake of simplicity, we consider an isothermal plasma. Figure R and Figure 



Table 1. Notations for compressible turbulence 



Notation 


Meaning 


a, c s , C/, Va 


sound, slow, fast, and Alfven speed 


5V, {SV) S , (5V) f , (SV)a 


random (rms) velocity 




Previously we used V for the rms velocity 


vi, {vi)s, (vi)f, (vi)a 


velocity at scale I 


Vk, (V k )s, (V k )/, (v k ) A 


velocity vector at wavevector k 


Bo (=ku), kL, k, §, ... 


unit vectors 


6,0 


displacement vectors 
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Fig. 8. Waves in real space. We show the directions of displacement vectors for a slow 
wave (left) and a fast wave (right). Note that £ s lies between d and Bo (= kn) and 
between k and kx- Again, 9 is perpendicular to k and parallel to the wave front. 
Note also that, for the fast wave, for example, density (inferred by the directions of the 
displacement vectors) becomes higher where field lines are closer, resulting in a strong 
restoring force, which is why fast waves are faster than slow waves. 



give schematics of slow and fast waves. For slow and fast waves, Bo, Vk (oc £), 
and k are in the same plane. On the other hand, for Alfven waves, the velocity 
of the fluid element (vu)a is orthogonal to the Bo — k plane. 

As before, the Alfven speed is Va = Bq/ where po is the average 
density. Fast and slow speeds are 



~ ja 2 + V% ± yj(a 2 + Vl) 2 - 4a 2 Vj cos 2 j ' . 



(25) 



where 6 is the angle between Bo and k. See Table [j] for the definition of other 
variables. When (3 ((3 — P g / Pb=2ci 2 /Vf ; P g — gas pressure, Pb— magnetic pres- 
sure; hereinafter f3 = average (3 = P 3 /Pb) goes to zero, we have 



c f « V A , 



acosO. 



(26) 



Figure |^ shows directions of displacement (or, directions of velocity) vectors 
for these three modes. We will call them the basis vectors for these modes. 
The Alfven basis is perpendicular to both k and B , and coincides with the 
azimuthal vector ^ in a spherical-polar coordinate system. Here hatted vectors 
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are unit vectors. The fast basis £/ lies between k and 
if « 



1-VD + /3/2 [k ' ' 



k 



fe||ki| + /c^kj 



(27) 



where D = (1 + /3/2) 2 - 2/3 cos 2 0, and /? is the averaged /3 (=P g /P B ). The slow 
basis £ s lies between and B (=ki|): 

1 - 

is oc fc||k|| 



1-VD- P/2\k 12 



1 + vD + (3/2 



ki ki 



(28) 



The two vectors £/ and £ s are mutually orthogonal. Proper normalizations are 
required for both bases to make them unit-length. 

When {3 goes to zero (i.e. the magnetically dominated regime), £/ becomes 
parallel to kj_ and £ s becomes parallel to Bo (Fig. 0b). The sine of the angle 
between Bo and £ s is (/3/2) sin0 cos 0. When (3 goes to infinity (i.e. gas pres- 
sure dominated regime )0, becomes parallel to k and £ s becomes parallel to 
(Fig. 0c). This is the incompressible limit. In this limit, the slow mode is 
sometimes called the pseudo- Alfven mode pa . 



6.2 Theoretical considerations 

Here we address the issue of mode coupling in a low (3 plasma. It is reasonable 
to suppose that in the limit where /5 > 1 turbulence for Mach numbers (M s = 
SV/a) less than unity should be largely similar to the exactly incompressible 
regime. Thus, Lithwick & Goldreich Q conjectured that the GS95 relations 
are applicable to slow and Alfven modes with the fast modes decoupled. They 
also mentioned that this relation can carry on for low (3 plasmas. For /3 > 1 
and M s > 1, we are in the regime of hydrodynamic compressible turbulence for 
which no theory exists, as far as we know. 

In the diffuse interstellar medium (3 is typically less than unity. In addition, it 
is ~ 0.1 or less for molecular clouds. There are a few simple arguments suggest- 
ing that MHD theory can be formulated in the regime where the Alfven Mach 
number (= SV/Va) is less than unity, although this is not a universally accepted 
assumption. Alfven modes describe incompressible motions. Arguments in GS95 
are suggestive that the coupling of Alfven to fast and slow modes will be weak. 
Consequently, we expect that in this regime the Alfven cascade should follow 
the GS95 scaling. Moreover the slow waves are likely to evolve passively Q . For 
a <^ Va their nonlinear evolution should be governed by Alfven modes so that 
we expect the GS95 scaling for them as well. The phase velocity of Alfven waves 
and slow waves depend on a factor of cos and this enables modulation of the 
slow waves by the Alfven ones. However, fast waves do not have this factor and 
therefore cannot be modulated by the changes of the magnetic field direction 

13 In this section, we assume that external mean field is strong (i.e. Va > {SV)) but 
finite, so that f3 —> oo means the gas pressure P g — > oo. 
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Fig. 9. Different regimes of MHD turbulence. We consider the compressible sub- 
Alfvenic regime (shaded region). In this figure, Va represents the total Alfven speed 
(= y/Bl + iB^I^/Wp). V = Vfiou, = SV. Cho & Vishniac (l9) argued that, even when 
the external field is weak, small scales can follow the GS95-like scaling in incompress- 
ible MHD regime. Similarly, for compressible case, the small scales of super-Alfvenic 
compressible turbulence are expected to fall in the sub-Alfvenic compressible regime. 
Moreover, winding of magnetic field by turbulence increases the magnetic field energy 
and the super-Alfvenic turbulence becomes more and more magnetically dominated 
with V/Va —> 1. Region B is the region where, within the density fluctuations, the 
velocities can get super-Alfvenic. In the figure, we used equation ( p3| ) to determine the 
borderline between region A and B. 



associated with Alfven waves. The coupling between the modes is through the 
modulation of the local Alfven velocity and therefore is weak. 

For Alfven Mach number (Ma) larger than unity a shock- type regime is 
expected. However, generation of magnetic field by turbulence |19 is expected 
for such a regime. It will make the steady state turbulence approach Ma ~ 1-F] 
Therefore in Cho & Lazarian [^l| we consider turbulence in the limit M s > 1, 
Ma < 1, and (3 < 1 (Fig. ^|). For these simulations, we mostly used M s ~ 2.2, 
Ma ~ 0.7, and (3 ~ 0.2. The Alfven speed of the mean external field is similar 
to the rms velocity (Va — 1,5V ~ 0.7, a = V0.1), and we used an isothermal 
equation of state. 

14 We suspect that simulations that show super-Alfvenic turbulence is widely spread 
in the ISM might not evolve for a long enough time to reach the steady state. 
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Fig. 10. (Left) Decay of Alfvenic turbulence. The generation of fast and slow waves is 
not efficient. (3 ~ 0.2, M s ~ 3. (ffipfti) The ratio of (SV)} to (<5V)a- The stronger the 
external field (Bo) is, the more suppressed the coupling is. The ratio is not sensitive to 
/3. From p| 



Although the scaling relations presented below are applicable to sub-Alfvenic 
turbulence, we cautiously speculate that small scales of super-Alfvenic turbu- 
lence might follow similar scalings. Boldyrev, Nordlund, & Padoan obtained 
energy spectra close to E(k) ~ fc -174 in solenoidally driven super-Alfvenic su- 
personic turbulence simulations. The spectra are close to the Kolmogorov spec- 
trum (~ fc~ 5 / 3 ), rather than shock-dominated spectrum (~ fc -2 ). This result 
might imply that small scales of super-Alfvenic MHD turbulence can be de- 
scribed by our sub-Alfvenic model presented below, which predicts Kolmogorov- 
type spectra for Alfven and slow modes. 



6.3 Coupling of MHD modes and Scaling of Alfven modes 



Alfven modes are not susceptible to collisionless damping (see 15611Cfl and ref- 



erences therein) that damps slow and fast modes. Therefore, we mainly consider 
the transfer of energy from Alfven waves to compressible MHD waves (i.e. to 
the slow and fast modes). 

In Cho & Lazarian plj , we carry out simulations to check the strength of the 
mode-mode coupling. We first obtain a data cube from a driven compressible 
numerical simulation with Bq/ ^/AiFpo — 1. Then, after turning off the driving 
force, we let the turbulence decay. We go through the following procedures be- 
fore we let the turbulence decay. We first remove slow and fast modes in Fourier 
space and retain only Alfven modes. We also change the value of Bo preserving 
its original direction. We use the same constant initial density po for all simula- 
tions. We assign a new constant initial gas pressure P g After doing all these 

15 The changes of both Bo and P g preserve the Alfven character of perturbations. 
In Fourier space, the mean magnetic field (Bo) is the amplitude of k = compo- 
nent. Alfven components in Fourier space are for k ^ and their directions are 
parallel/anti-parallel to £a (= Bo x kj_). The direction of £a does not depend on the 
magnitude of So or P g . 
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Fig. 11. (a) Alfven spectra follow a Kolmogorov-like power law. (b) The second-order 
structure function (SF2 =< v(x + r) — v(x) >) for Alfven velocity shows anisotropy 
similar to the GS95. Conturs represent eddy shapes. From H]. 



procedures, we let the turbulence decay. We repeat the above procedures for 
different values of Bq and P g . Fig. [To| a shows the evolution of the kinetic energy 
of a simulation. The solid line represents the kinetic energy of Alfven modes. It 
is clear that Alfven waves are poorly coupled to the compressible modes, and 
do not generate them efficiently |^ Therefore, we expect that Alfven modes will 
follow the same scaling relation as in the incompressible case. Fig. p~0| b shows 
that the coupling gets weaker as Bq increases: 



(SV)} (6VU 

mi* b ■ 



(29) 



The ratio of (6V) 2 S to (SV)\ is proportional to (5V)\/B$. 

This marginal coupling is in good agreement with a claim in GS95, as well as 
earlier numerical studies where the velocity was decomposed into a compressible 
component vc and a solenoidal component vs. The compressible component is 
curl-free and parallel to the wave vector k in Fourier space. The solenoidal com- 
ponent is divergence-free and perpendicular to k. The ratio \ = (SV)c/{SV)s 
is an important parameter that determines the strength of any shock 133| 141] 



Porter, Woodward, & Pouquet [ 138 1 performed a hydrodynamic simulation of 
decaying turbulence with an initial sonic Mach number of unity and found that 
X 2 evolves toward ~ 0.11. Matthaeus et al. ]105|| carried out simulations of de- 



caying weakly compressible MHD turbulence [ 183 1 and found that \ 2 ~ 0(M 2 ), 
where M s is the sonic Mach number. In [|l0| a weak generation of compressible 
components in solenoidally driven super- Alfvenic supersonic turbulence simula- 
tions was obtained. 

Fig. ^ shows that the spectrum and the anisotropy of Alfven waves in this 
limit are compatible with the GS95 model: 

Spectrum of Alfven Modes: E(k) on k ± 5 ^ 3 , (30) 

16 As correctly pointed out by Zweibel (this volume) there is always residual coupling 
between Alfven and compressible modes due to steepening of Alfven modes. However, 
this steepening happens on time-scales much longer than the cascading time-scale. 
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Fig. 12. (a) Slow spectra also follow a Kolmogorov-like power law. (b) Slow modes 
show anisotropy similar to the GS95 theory. From km. 
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and scale-dependent anisotropy fc|| cx fcj_' that is compatible with the GS95 
theory. 

6.4 Scaling of the slow modes 

Slow waves are somewhat similar to pseudo-Alfven waves (in the incompressible 
limit). First, the directions of displacement (i.e. £ s ) of both waves are similar 
when anisotropy is present. The vector £ s is always between 9 and k|| . In Figure 
, we can see that the angle between 9 and k|| gets smaller when fc|| <C k±. 
Therefore, when there is anisotropy (i.e. fc|| <C k±), £ s of a low (3 plasma be- 
comes similar to that of a high /3 plasma. Second, the angular dependence in the 
dispersion relation c s ~ a cos 9 is identical to that of pseudo-Alfven waves (the 
only difference is that, in slow waves, the sound speed a is present instead of the 
Alfven speed Va)- 

Goldreich & Sridhar |^(| argued that the pseudo-Alfven waves are slaved to 
the shear- Alfven (i.e. ordinary Alfven) waves in the presence of a strong B , 
meaning that the energy cascade of pseudo-Alfven modes is primarily mediated 
by the shear- Alfveen waves. This is because pseudo-Alfven waves do not provide 
efficient shearing motions. Similar arguments are applicable to slow waves in a 
low (3 plasma (see also for high-/? plasmas). As a result, we conjecture 
that slow modes follow a scaling similar to the GS95 model : 

Spectrum of Slow Modes: E s (k) oc k ± 5 ^ 3 . (31) 

Fig. |l^ a shows the spectra of slow modes. For velocity, the slope is close 
to —5/3. Fig. [12) 3 shows the contours of equal second-order structure function 
(SF2) of slow velocity, which are compatible with fc|| cx k 2 ^ 3 scaling. 

In low P plasmas, density fluctuations are dominated by slow waves plfl . 
From the continuity equation p — p\7 ■ v 



ujp k = p a k ■ v fe , 



(32) 
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we have, for slow modes, (pk)s ~ Po(vk)s/o- Hence, this simple argument implies 

'Sp\ _ (5V) S 



P 



M s , (33) 



where we assume that (5V) S ~ (SV)a and M s is the Mach number. On the 
other hand, only a small amount of magnetic field is produced by the slow 
waves. Similarly, using the induction equation (wbfc = kx (Bo X v*,)), we have 

which means that equipartition between kinetic and magnetic energy is not guar- 
anteed in low j3 plasmas. In fact, in Fig. |l2|a, the power spectrum for density 
fluctuations has a much larger amplitude than the magnetic field power spec- 
trum. Since density fluctuations are caused mostly by the slow waves and mag- 
netic fluctuation is caused mostly by Alfven and fast modes, we do not expect 
a strong correlation between density and magnetic field, which agrees with the 



ISM simulations [130 126,170 1 



6.5 Scaling of the fast modes 

Figure |l3| shows fast modes are isotropic. The resonance conditions for interact- 
ing fast waves are: 

ui + UJ2 = w 3 , (35) 
k!+k 2 =k 3 . (36) 

Since uj oc k for the fast modes, the resonance conditions can be met only when 
all three k vectors are collinear. This means that the direction of energy cascade 
is radial in Fourier space, and we expect an isotropic distribution of energy in 
Fourier space. 

Using the constancy of energy cascade and uncertainty principle, we can 
derive an IK-like energy spectrum for fast waves. The constancy of cascade rate 
reads 

— l — = - = constant. (37) 

teas teas 

On the other hand, t cas can be estimated as 

t „ v * V k f, sl 

( V ' Vv) k Ep+q=k kv P V q 

If contributions are random, the denominator can be written by the square root 
of the number of interactions (VN) times strength of individual interactions 
(~ kv%) Here we assume locality of interactions: p ~ q ~ k. Due to the 

17 To be exact, the strength of individual interactions is ~ fci;^sin6l, where 6 is the 
angle between k and Bo. Thus marginal anisotropy is expected. It will be investigated 
elsewhere. 
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uncertainly principle, the number of interactions becomes Af ~ k(Ak) 2 , where 
Ak is the typical transversal (i.e. not radial) separation between two wave vec- 



tors p and q (with p + q = k). Therefore, the denominator of equation (38) is 
(k(Ak) 2 ) 1 ^ 2 kv 2 . We obtain an independent expression for t cas from the uncer- 
tainty principle (t cas Auj ~ 1 with Au> ~ Ak{Ak/k)). From this and equation 
@, we get t cas ~ t 1 J a 2 s /{k 2 v k ) 1 which yields 

t cas ~ l/k 4 V 2 k . (39) 



Combining equations (37) and ([39]), we obtain v 2 ^ k 7 / 2 , or E? (k) ~ fc 2 ^ 2 . 



fc 3 ' 2 . This is very similar to acoustic turbulence, turbulence caused by inter- 



acting sound waves 

pqf8ip 7|]. Zakharov & Sagdeev |18| found E(k) oc k~ z l 2 . 
However, there is debate about the exact scaling of acoustic turbulence. Here we 
cautiously claim that our numerical results are compatible with the Zakharov & 
Sagdeev scaling: 

Spectrum of Fast Modes: E f (k) - k~ 3/2 . (40) 

Magnetic field perturbations are mostly affected by fast modes |JlJ when (3 
is small: 

if (SV) A - (SV) S . 

The turbulent cascade of fast modes is expected to be slow and in the absence 
of collisionlcss damping they are expected to propagate in turbulent media over 
distances considerably larger than Alfven or slow modes. This effect is difficult to 
observe in numerical simulations with AB ~ B . A modification of the spectrum 



in the presence of the collisionless damping is presented in [176 
7 Astrophysical Implications 



Many astrophysical problems require some knowledge of the scaling properties of 
turbulence. Therefore we expect a wide range of applications of the established 
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scaling relations. Here we show how recent breakthroughs in understanding MHD 
turbulence affect a few selected issues. 



7.1 Cosmic ray propagation 

The propagation of cosmic rays is mainly determined by their interactions with 
electromagnetic fluctuations in interstellar medium. The resonant interaction of 
cosmic ray particles with MHD turbulence has been repeatedly suggested as the 
main mechanism for scattering and isotropizing cosmic rays. In these analysis it 
is us ually assumed that the turbulence is isotropic with a Kolmogorov spectrum 



(see [151 ). How should these calculations be modified? 

The essence of the mechanism is rather simple. Particles moving with velocity 
v interact with a resonant Alfven wave of frequency u) = k\\Vji-\-nQ (n = ±1, 2...), 
where Q — Qq/j is the gyrofrequency of relativistic particles, /i is the cosine of 
the pitch angle. From the resonant condition above, we know that the most 
important interaction occurs at fcii ~ Q/vfj, ~ (/nr^) , where is Larmor 
radius of the high-ener gy p articles. 



The calculations in [ 178 1 that made use of tensor (|l4|) provided the scattering 



efficiency of anisotropic Alfvenic turbulence. The results are compared in Fig. ft4| a 
with the predictions of the scattering on isotropic Kolmogorov-type magnetic 
fluctuations and also with earlier calculations by Chandran Jl7j . The latter used 
rather ad hoc form of the tensor to describe magnetic fluctuations within the 
Goldreich-Sridhar theory. We see from Fig |l4| a that the scattering is substantially 
suppressed, compared to the Kolmogorov turbulence that is usually used for 
scattering calculations. This happens, first of all, because most turbulent energy 
in GS95 turbulence goes to k± so that there is much less energy left in the 
resonance point fcy = (/w^,) -1 . Furthermore, k±_ 3> kn means k± 3> r^ 1 so that 
cosmic ray particles cover many eddies during one gyration. This random walk 
decreases the scattering efficiency by a factor of (f2/k±v±)z = {tl/Ix)^, where 
l± is the turbulence scale perpendicular to magnetic field. 

Thus the gyroresonance with Alfvenic turbulence is not an effective scattering 
mechanism for cosmic rays if turbulence is injected on the large scales, since the 
degree of anisotropy increases on smaller scales. However, if energy is injected 
isotropically at small scales, the resulting turbulence would be more isotropic 
and scattering will be more efficient. Scattering by undamped fast modes is more 



efficient than the Kolmogorov theory would predict. Yan & Lazarian | 178 | per- 
formed calculations taking into account the collisionless damping of fast modes 
and showed that the gyroresonance scattering by fast modes is the dominant 
scattering mechanism. 

There is another important property of turbulence that was neglected in 
earlier work. When cosmic rays stream at a velocity much larger than Alfven 
velocity, they can excite resonant MHD waves, which in turn scatter cosmic rays. 
This is the 'streaming instability'. It is usually assumed that this instability can 
confine cosmic rays with energies less than lOOGeV However, this is true 
only in an idealized situation when there is no background MHD turbulence. As 
noted earlier, the rates of turbulent decay are very fast and excited perturbations 
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Grain velocities induced by turbulence 




Fig. 14. Applications. (By Lazarian & Yan). Left: Cosmic ray scattering for realistic 
MHD turbulence is reduced substantially compared to scattering by isotropic turbu- 
lence, but still larger than estimates in Chandran Il7j. Right: The dust acceleration by 
turbulence is reduced compared to the accepted estimates in Draine p3[. 



should vanish quickly. In | 178 | we find that the streaming instability is only 
applicable to particles with energies < O.lUGeV, which is less than the energy 
of most cosmic ray particles. This result casts doubt on the self-confinement 
mechanism discussed by previous authors. 

All these findings tend to support the alternative picture of cosmic ray diffu- 
sion advocated by Jokipii (see |69|| ). In this picture cosmic rays follow magnetic 
field lines, but the magnetic field wanders. The rate of this wandering can be 
calculated from the established turbulence scaling laws. 

Knowledge of the scattering rates is essential for understanding both the first 
order and the second order Fermi acceleration. The first order Fermi acceleration, 
may be important for a wide range of phenomena from clusters of galaxies and 
gamma-ray bursts (GRBs) to solar flares. Results obtained in [178 where the 
discovered properties of MHD turbulence are used proved to be very different 
from earlier estimates. 



7.2 Grain dynamics 

Turbulence induces relative dust grain motions and leads to grain-grain colli- 
sions. These collisions determine grain size distribution, which affects most dust 
properties, including absorption and H2 formation. Unfortunately as in the case 
of cosmic rays, earlier work appealed to hydrodynamic turbulence to predict 
grain relative velocities (see |74 

The differences between the hydrodynamic and MHD calculations stem from 
(a) grain charges, which couple grains to the magnetic field, (b) the anisotropy 
of MHD cascade, and (c) the direct interaction of charged grains with magnetic 
perturbations. Effects (a) and (b) are considered in Lazarian & Yan |92|| , while 
(c) is considered in Yan & Lazarian (2002b; in preparation). As consequence the 
picture of grain dynamics is substantially altered. 



EM, 33, 124, 167, 3£ 
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Consider grain charge first. If a grain's Larmor period tl — 2 / nm gr c/ qB is 
shorter than the gas drag time tdrag, grain perpendicular motions are constrained 
by magnetic field. Their velocity dispersion is determined by the turbulence eddy 
whose turnover period is ~ tl instead of the drag time Q . 

Accounting for the anisotropy of MHD turbulence it is convenient to consider 
separately grain motions parallel and perpendicular to magnetic field. The per- 
pendicular motion is influenced by the Alfven modes, which have a Kolmogorov 
spectrum. The parallel motion is subjected to compressible modes which scale 

— 1/2 

as v\\ oc ku . In addition we should account for viscous forces. When the eddy 

turnover time is of the order of tdamp ~ v~ 1 kJ 2 , the turbulence is viscously 
damped. Thus grains sample only a part of the eddy before gaining the velocity 
of the ambient gas if tx or tdrag < tdamp- The results are shown in Fig. |l4|b. 

The direct interaction of the charged grains with turbulent magnetic field 
results in a stochastic acceleration that can potentially provide grains with su- 
personic velocities. 



7.3 Turbulence in HII regions 

Lithwick & Goldreich |)6| addressed the issue of the origin of density fluctu- 
ations within HII regions. There the gas pressure is larger than the magnetic 
pressure (the 'high beta' regime) and they conjectured that fast waves, which 
are essentially sound waves, would be decoupled from the rest of the cascade. 
They found that density fluctuations are due to the slow mode and the entropy 
mode, which are passively mixed by shear Alfven waves and follow a Kolmogorov 
spectrum. They also found that slow mode density fluctuations are proportional 
to l/vT?- On the other hand, the entropy mode density fluctuations are sup- 
pressed when cooling is faster than the cascade time. Lithwick & Goldreich 
also gave detailed discussions about density fluctuations on various scales in the 
ISM, e.g. proton gyro-radius. These results are important as radio-wave scin- 
tillation observations can constrain the nature of MHD turbulence in the ISM, 
especially in the HII regions. Lithwick & Goldreich argued that the turbu- 
lent cascade survives ion-neutral damping only when a high degree of ionization 
is present. However, the study by Cho, Lazarian & Vishniac [Q suggests that 
the magnetic fluctuations protrude below the damping scale and the results of 
[ p6| should be revised. 



7.4 Tiny-Scale Atomic Structures 

The intermittent small scale structures in § |5.2| should have important implica- 
tions for transport processes (heat, cosmic rays, etc.) in partially ionized plasmas. 
We also speculate that they might have some relation to the tiny-scale atomic 
structures (TSAS). Heiles introduced the term TSAS for the mysterious 
H I absorbing structures on scales from thousands to tens of AU, discovered by 
Dieter, Welch & Romney J3i|, Analogs are observed in Nal and Call |l08p| [j] 

[28| | an- 



and in molecular gas [ 102 1 . Recently Deshpande, Dwarakanath & Goss 
alyzed channel maps of opacity fluctuations toward Cas A and Cygnus A. They 
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found that the amplitudes of density fluctuations at scales less than 0.1 pc are 
far larger than expected from extrapolation from larger scales, possibly explain- 
ing TSAS. This study, however, cannot answer what confines those presumably 
overpressured (but very quiescent!) blobs of gas. Deshpande |27j related those 
structures to the shallow spectrum of interstellar turbulence. 

Figure [|b indicates that while velocity decreases rapidly, but not exponen- 
tially, below the viscous damping scale, the magnetic field fluctuations persist, 
thereby providing nonthermal pressure support. Magnetic structures perpendic- 
ular to the mean magnetic field are compensated by pressure gradients. Our 
calculations so far are produced using incompressible code In the case of 
compressible media, wc expect the pressure fluctuations to entail density fluctu- 
ations reminiscent of the Deshpande et al. [^8| observations. 

The calculations in Cho, Lazarian & Vishniac are applicable on scales 
from the viscous damping scale (determined by equating the energy transfer 
rate with the viscous damping rate; ~ 0.1 pc in the Warm Neutral Medium 
with n = 0.4 cm~ 3 , T= 6000 K) to the ion- neutral decoupling scale (the scale 
at which viscous drag on ions becomes comparable to the neutral drag; <C 0.1 
pc). Below the viscous scale the fluctuations of magnetic field obey the damped 
regime shown in Figure ||b and produce density fluctuations. For typical Cold 
Neutral Medium gas, the scale of neutral-ion decoupling decreases to ~ 70AU, 
and is less for denser gas. TSAS may be created by strongly nonlinear MHD 
turbulence! 

A simple technique of estimating magnetic field was suggested by Chan- 
drasekhar & Fermi |l8j (see also review by Ostriker, this volume). According 
to it, the fluctuations of magnetic field that can be measured from polarization 
maps are related to velocity fluctuations measured through Doppler broadening 
Sb/y/Anp ~ 5v. The existence of the damped regime of MHD turbulence suggests 
that this technique is not applicable to very small scales in partially ionized gas. 

7.5 Magnetic Reconnection 

Magnetic reconnection is the fundamental process that allows magnetic fields 
to change their topology, despite being 'frozen' into highly conducting plasmas. 
It is the key process for solar flares, the magnetic dynamo, the acceleration of 
energetic particles, etc. According to the Lazarian & Vishniac model |8!| (see 
also review |H)) and Vishniac, Lazarian, & Cho, this volume) of stochastic recon- 
nection, this process is controlled by the turbulent wandering of magnetic field. 
The exact properties of the turbulent cascade are especially important for the 
viscously damped regime present in partially ionized gas. However, it is shown in 
Lazarian, Vishniac & Cho |n) that the reconnection rates are sufficiently high 
in this case. The implications of the finding for the removal of magnetic flux 
during star formation is to be evaluted yet. 
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7.6 Support and Compression of Molecular Clouds 

To understand the dynamics of molecular clouds and star formation it is nec- 



essary to understand turbulence. In a recent review [ 106 1 McKee pointed out 
that the fast damping of MHD turbulence observed in numerical simulations is 
difficult to reconcile with the fact that "a GMC such as G216, which has no 
visible star formation, can have a level of turbulence that exceeds that in the 
Rosette molecular cloud, which has an embedded OB association" . He pointed 
out that the conclusions obtained on the basis of numerics should be treated 
with caution as they do not resolve the microscales. 

In typical astrophysical conditions the sources of turbulence are localized 
both in space (stars) and time (stellar outflows; supernovae), and the outgoing 
waves have much larger amplitudes than the background waves (we call this 
situation "imbalanced cascade" ). Fig. ||a shows that the turbulent damping could 
be substantially reduced in this situation. Moreover, even in a balanced regime, 
we expect fast modes to be subjected to slow non-linear damping. 

At the same time, it worth mentioning, that turbulence can not only sup- 
port, but also compress molecular clouds. Clouds can be compressed by external 
turbulence feeding into them and depositing energy and momentum. Myers & 



Lazarian [12C] explained observed infalling motions of molecular gas surround- 



ing dense cores [164,94] in this way, based on ion-neutral damping. The infall 
rate is proportional to the rate of turbulence damping. Therefore, fast non- linear 
damping associated with the Alfvenic turbulence should enhance the infall. 



7.7 Heating of Diffuse Ionized Gas 



The "Diffuse Ionized Gas" (DIG), or equivalently the "Reynolds layer" within 
the Milky Way, is detected by rather faint but ubiquitous Galactic Ha emission 
[144 145]. Such emission is found in several other spirals as well [142,143 127,128]. 
In the Galaxy, the Reynolds layer contains a substantial portion of the H + in 
the ISM. Current models generally involve photoionization from the OB stars, 
although how the Lyman continuum radiation from OB stars can penetrate the 
neutral H layer remains controversial. 

The observations show strong [S II] A6717 and [N II] A6583 that increase rel- 
ative to Ha with distance z above the planes of various galaxies, including the 
Milky Way 145]. The only reasonable conclusion is that there is an additional 
source of heating in the ISM that dominates over photoionization heating at low 
densities. It has been proposed that carbonaceous molecules provide the excess 
heating through the photoelectric effect [177|, but this explanation is not unique. 
Heating by turbulence, surely present, may dominate. Minter & Spangler flllCfl 
suggested a heating rate that is adequate to explain the [S II] /Ha and [N II] /Ha 
ratios, but did not take nonlinear interactions into account, thereby underesti- 
mating the heating. A new study that would capitalize on the new understanding 
of MHD turbulence (damping of Alfven and fast modes, imbalance etc.) is on 
our agenda. 
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8 Observational Tests 



Comparing numerics with observations is a challenging problem. Ostriker (this 
volume) discusses PDF, clump identification, and linewidth-scale relations as 
possible diagnostics and outlines problems with any of those approaches. A use 
of spectral line data cubes and application to it of different techniques (e.g. spec- 
tral correlation function, principal component analysis, wavelets, etc.) can be 
found in the review [f79| . Here we shall concentrate on comparing spectra from 
observations with our theoretical expectations (see also review J86|). 



8.1 Is the Big Power Law real? 

We have mentioned above that observations suggest that the Kolmogorov power 
law should span from AU to kpc scales. Kolmogorov scaling is exactly what one 
would expect from the GS95 picture when the observations sample magnetic field 
in the system of reference aligned with the mean magnetic field (see Fig. [l^ for 
examples of observable quantities). Indeed, it is obvious from Fig. |l| in the this 
system of reference (i.e. global system of reference) the locally defined scalings 
of fc|| with k± are not valid. Indeed, one can easily see that the fluctuations 
perpendicular to the local direction of magnetic field dominate both the statistics 
measured perpendicular to the mean magnetic field and parallel to it. As the 
result in the reference system aligned with the mean magnetic field k'^ cx k' ± (see 

Fig. |l5|a,b,c) and according to equation ( pi] ) E{k) ~ fc -5 / 3 will be measured. 

Ambiguities in measurements reviewed in 0], however, make it uncertain 
whether or not the Big Power Law should be taken at face value. Still, we note 
that agreement between an observed power spectrum and theoretical expecta- 
tions is far more significant than just an approximate fit between observations 
and numerics. The latter is definitely not unique and is a priori suspect in view 
of the huge difference in terms of Re and Rm between any numerical simulation 
and the ISM. 

To test the Big Power Law properly, it is important to extend the theory 
of Velocity Channel Analysis (VCA) Q,^6| by including self-absorption in the 
analysis of turbulent emission lines. It can then be applied to regions of HI in 
the inner Galaxy Q or CO ^63^16^ 1^. Lazarian & Po gosyan Q provided 



results consistent with observations. For instance, the study predicts that in 
the presence of absorption the emitted power in the line is proportional to k~ 3 , 
exactly what is seen in Dickey et al. |30| for HI in the inner region of the Galactic 
plane. 

An application of the VCA to different emission lines (e.g. Ha, [N II], [S II] 
etc.) would help to answer the question of whether or not ISM turbulence is 

18 Brunt & Heyer jjlj applied Principal Component Analysis (see (5^]), to simulated 
CO data and found empirical relations between the statistics of velocities, eigen- 
vectors and eigenimages. However, they note that somehow their relation does not 
depend on the absorption coefficient for the limited range of absorptions they tested. 
Unfortunately, their analysis does not seem to be applicable to Galactic HI and the 
applicability of their technique to the correlated velocity and density fields is unclear. 
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a large scale cascade with various phases of the ISM interconnected through a 
dynamically important magnetic field. A contrasting possibility is that various 
phases form their own cascades. 

HI is rather smoothly distributed across the sky. Therefore, the effects of 
image edges does not pose a problem for the statistical analysis. Results by 
Stanimirovic (private communication) show that the Fourier analysis of the SMC 
image and a more laborious wavelet analysis]^ provide identical results. However, 
when dealing with molecular clouds we might expect that t he c loud edges become 



important. Therefore the incorporation of wavelets (see [184 ) within the VCA 



is a natural step to make. The only difference would be to apply the wavelets 
instead of Fourier transforms to the channel maps. 

Synchrotron fluctuations and fluctuations of polarized radiation arising from 
aligned dust should be used to study magnetic field statistics. As we mentioned 
earlier Cho & Lazarian J2^| has shown that those fluctuations are consistent with 
Kolmogorov scaling. More studies in this direction are necessary. The fact that 
those fluctuctuations interfere with the CMB studies garantees that in the near 
future we shall have a lot of relevant data. 

One should remember, however, that the measured power index of fluctua- 
tions may not correspond to the spectral index of the underlying turbulence. For 
instance, it is shown in |159| that while the actual turbulence in SMC is close 
to being Kolmogorov, depending on the thickness of the slice, the spectral in- 
dex of intensity fluctuations within channel maps span from ~ —2.8 to ~ —3.4. 
Similarly, it is shown in p2| that for Kolmogorov turbulence the spectral in- 
dex of observable fluctuations may vary from ~ —1 to —3.7 depending on how 
observations sample turbulence. 



8.2 Does turbulence reveal magnetic field direction? 

Anisotropy of Alfvenic turbulence is a definite prediction of the GS95 theory. We 
mentioned earlier that the scale-dependent anisotropy can only be revealed in the 
local frame of reference, which in practical terms require direct measurements, 
e.g. with spacecraft. Measurements of the Solar wind magnetic turbulence have 
failed so far to reveal the differential scaling of the turbulence in terms of feii and 
k± , but these measurements are inconclusive |37| . If measurement are performed 
in a global system of reference, as is the case with observations, they should reveal 
anisotropy in the direction of the mean magnetic field. 

In isotropic turbulence, correlations depend only on the distance between 
sampling points. Contours of equal correlation are circular in this case. The 
presence of a magnetic field introduces anisotropy and these contours become 
elongated with a symmetry axis given by the magnetic field. To study turbulence 
anisotropy, we can measure contours of equal correlation corresponding to the 
data within various velocity bins. The results obtained with simulated data are 
shown in Fig. |l5|. 

19 Wavelet analysis involves determining the deviation of each pixel from a weighted 
average of the pixels at a particular projected distance from it. 
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Fig. 15. Observational test of the synthetic data by Cho, Esquivel & Lazarian. Con- 
tours of equal correlation obtained with Centroids of Velocity ( (a)) and with Spectral 
Correlation Function (SCF)((b)). The direction of anisotropy reveals the direction of 
projected magnetic field. Combined with the anisotropy analysis, the SCF (introduced 
by Alyssa Goodman) is likely to become even more useful tool, (c) Contours of equal 
correlation obtained for synthetic synchrotron intensity map. ( d) shows the 2D genus of 
the Gaussian distribution (smooth analytical curve) against the genus for the isother- 
mal compressible MHD simulations with Mach number ~2.5 (dotted curve). 



Since the degree of anisotropy is related to the strength of the magnetic 
field, studies of anisotropy can provide the means to analyze magnetic fields. It 
is important to study different data sets and channel maps for the anisotropy. 
Optical and infrared polarimctry can benchmark the anisotropies in correlation 
functions. We hope that anisotropies will reveal magnetic field structure within 
dark clouds where grain alignment and therefore polarimetry fails (see for a 
review of grain alignment). 

Not only velocity statistics can be used for such an analysis. Lazarian & 
Shutenkov (see review |76j) showed that the mean magnetic field must lead 
to anisotropies in the synchrotron statistics. Lazarian & Chibisov |8l| pointed 
out that using HI regions as screens for radiation at the decameter wavelength 
it should be possible to study the 3D distribution of the magnetic field. Fig. |l5| c 
shows the anisotropy of synchrotron statistics available through simulations. 
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8.3 How else can we compare observations and simulations? 

Velocity and density power spectra do not provide a complete description of 
turbulence. Intermittency (variations in the strength of the turbulent cascade) 
and its topology in the presence of different phases are not described by the 
power spectrum. Use of the higher moments is possible (see discussion of the 3 
point statistics, or bispectrum, in Lazarian |79||), but is limited by the noise in 
the observational data. 

"Genus analysis" is a good tool for studying the topology of turbulence (see 
the review |8^| ). This tool has already been successfully applied to cosmology 
[ jjil) . Consider an area on the sky with contours of projected density. The 2D 
genus, G(^), is the difference between the number of regions with a projected 
density higher than v and those with densities lower than v. Fig. p"5| d shows the 
2D genus as the function of v for a Gaussian distribution of densities (completely 
symmetric curve), for MHD isothermal simulations with Mach number ^2.5. It 
is shown in [^6| that the genus of the Small Magellanic Cloud is very different 
from that in Fig. [l5|d, while the spectra in both cases are similar. 

9 Summary 

Recently there have been significant advances in the field of MHD turbulence: 

1. The first self-consistent model (GS95) of incompressible MHD turbulence 
that is supported by both numerical simulations and observations is now 
available. The major predictions of the model are scale-dependent anisotropy 

tx fc^ 3 ) and a Kolmogorov energy spectrum (E(k) oc fc -5 / 3 ). 

2. There have been substantial advances in understanding compressible MHD. 
Simulations of compressible MHD turbulence show that there is a weak cou- 
pling between Alfven waves and compressible MHD waves and that the 
Alfven modes follow the Goldreich-Sridhar scaling. Fast modes, however, 
decouple and exhibit isotropy. 

3. Contrary to general belief, in typical interstellar conditions, magnetic fields 
can have rich structures below the scale at which motions are damped by the 
viscosity created by neutral drag (the ambipolar diffusion damping scale). 

4. These advances will have a dramatic impact on our understanding of many 
fundamental interstellar processes, like cosmic-ray propagation, grain dy- 
namics, turbulent heating and molecular cloud stability. 

5. New techniques, e.g. VCA, allow observational tests of the theory. 
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